AFRL-SN-RS-TR-2004-44 
Final  Technical  Report 
February  2004 


STUDY  OF  LINEARIZATION  OF  OPTICAL 
POLYMER  MODULATORS 


University  of  Minnesota 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


AIR  FORCE  RESEARCH  LABORATORY 
SENSORS  DIRECTORATE 
ROME  RESEARCH  SITE 
ROME,  NEW  YORK 


STINFO  FINAL  REPORT 


This  report  has  been  reviewed  by  the  Air  Foree  Researeh  Laboratory,  Information 
Direetorate,  Publie  Affairs  Offiee  (IFOIPA)  and  is  releasable  to  the  National  Teehnieal 
Information  Serviee  (NTIS).  At  NTIS  it  will  be  releasable  to  the  general  publie, 
ineluding  foreign  nations. 


AFRL-SN-RS-TR-2004-44  has  been  reviewed  and  is  approved  for  publieation. 


APPROVED;  /s/ 

BRIAN  F.  MCKEON 
Project  Engineer 


EOR  THE  DIRECTOR:  /s/ 

RICHARD  G.  SHAUGHNESSY,  Et  Col,  USAE 
Chief,  Rome  Operations  Office 
Sensors  Directorate 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  074-0188 


Public  reporting  burden  for  this  coilection  of  information  is  estimated  to  average  1  hour  per  response,  inciuding  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  compieting  and  reviewing  this  coilection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  coliection  of  information,  inciuding 
suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Ariington,  VA  22202-4302, 
and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503 _ 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 


FEBRUARY  2004 


FINAL 


4.  TITLE  AND  SUBTITLE 


STUDY  OF  LINEARIZATION  OF  OPTICAL  POLYMER  MODULATORS 


6.  AUTHOR(S) 

Anand  Gopinath,  Ross  Schermer,  Jaesang  Oh,  Desalegn  Bereka 


Dec  02  -  Sep  03 


5.  FUNDING  NUMBERS 

G  -F30602-03-1-0013 
PE  -63175C 
PR  -6015 
TA  -  SN 
WU  -  01 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Minnesota 

University  of  Minnesota  Sponsored  Projects 

200  Oak  Street  SE 

Minneapoiis  MN  55455-2070 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFRL/SNDP 
26  Electronic  Parkway 
Rome  NY  13441-4514 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 

AFRL-SN-RS-TR-2004-44 


11.  SUPPLEMENTARY  NOTES 

AFRL  Project  Engineer:  Brian  F.  McKeon/SNDP/(315)  330-7348 

Brian.McKeon@rl.af.mil 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 

13.  ABSTRACT  (Maximum  200  Words) 

To  improve  the  Spur  Free  Dynamic  Range  of  analog  electro-optic  modulators  in  the  >10  GHz  regime,  techniques  for 
improving  the  linearity  of  these  devices  must  be  developed.  This  report  discusses  an  investigation  into  electro-optic 
directional  couplers  that  use  variable  coupling  in  polymer-based  waveguide  structures.  The  approach  begins  with 
determining  the  desired  response  function,  in  this  case  a  trapezoid.  The  Fourier  transform  method  is  then  applied  to 
synthesize  the  variable  coupling  function  which  results  in  the  desired  output  response.  This  initial  coupling  function  is 
iteratively  modified  using  Newton’s  method  to  obtain  a  coupling  design  that  can  be  implemented  in  a  polymer 
waveguide.  Approaches  to  achieving  the  variable  coupling,  including  regions  of  negative  coupling,  are  discussed. 


14.  SUBJECT  TERMS 

Electro-Optic  Modulator,  Optical  Polymer  Modulators,  Directional  Coupler,  Polymer  Optical 
Waveguide,  Mach-Zehnder,  Linearization  Techniques,  Spur  Free  Dynamic  Range 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


UNCLASSIFIED 


NSN  7540-01-280-5500 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18 
298-102 


Table  of  Contents 


1.  Objective  1 

2.  Introduction  1 

3.  Design  Technique  for  Coupler  Modulators  1 

4.  Calculation  of  Switching  Voltage-Length  Products  for  different  modulator  designs  5 

5.  Distortion  Calculations  12 

6.  Linearized  Coupler  Designs  in  Polymer  Waveguides  14 

6.1  Realization  of  Variable  Coupling  Designs  in  Polymer  Waveguides  15 

6.1.1  Variable  Coupling  Schemes  15 

6.1.2  Negative  Coupling  by  Phase  Shifts  1 7 

6.1.3  Electrode  Structure  for  Polymer  Coupler  17 

7.  Conclusions  18 

8.  References  19 


1 


List  of  Figures 


1 .  Schematic  of  the  uniform  coupling  coupler  modulator. 

2.  The  desired  trapezoidal  response  function  obtained  after  several  iterations,  the  normalized 
intensity  response  plotted  against  the  normalized  switching  voltage. 

3.  The  iterated  coupling  function  for  the  response  function  in  figure  2. 

4.  Plot  of  normalized  power  against  normalized  distance  along  the  coupler  for  different  bias 
values.  The  power  in  the  second  guide  of  the  coupler  with  trapezoidal  response  shown  in 
figure  2.  The  blue  curve  is  for  zero  bias,  green  is  bias  which  obtains  half  power,  and  the 
red  curve  is  bias  when  no  transfer  takes  place  into  the  second  guide. 

5.  Schematic  diagram  of  the  Mach  Zehnder  modulator. 

6.  The  Mach  Zehnder  calculated  response. 

7.  Schematic  diagram  of  the  uniform  coupler. 

8.  The  response  of  the  uniform  coupler  in  which  the  normalized  transmitted  power  is  plotted 

against  the  normalized  voltage.  The  value  of  p  for  this  coupler  is  VS  /2  for  the  push-pull 
case. 

9.  The  linear  response  coupler  modulator  coupling  function  obtained  from  the  Fourier 
transform  of  the  desired  linear  response. 

10.  The  linear  response  coupler  modulator  showing  the  response  function  as  obtained  from 
the  coupling  function  in  figure  9. 

1 1 .  The  trapezoidal  response  coupler  design  showing  the  coupling  function. 

12.  The  trapezoidal  response  coupler  design  showing  the  response  function,  desired  and 
calculated. 

13.  The  trapezoidal  response  coupler  design  from  figure  1 1  showing  the  modified  coupling 
function  to  have  flat  top  regions. 

14.  The  trapezoidal  response  coupler  design  from  figure  13  showing  the  resultant  response 
function. 

15.  The  constant  |k|  coupler  modulator  design  but  with  sign  changes  for  the  coupling. 

16.  The  constant  |k|  coupler  modulator  but  with  sign  changes  for  the  coupling.  The  response 
function  is  shown,  the  value  of  p  is  1.8  and  the  L/Lc ratio  is  2.0. 

1 7.  The  response  function  of  the  constant  |k|  coupler  modulator  design  but  with  sign  changes 
for  5  along  the  length  of  the  coupler. 

18.  The  response  function  of  the  constant  |k|  coupler  modulator  but  with  sign  changes  for 
5  along  the  length  of  the  coupler,  the  value  of  p  is  2.7  and  L/Lc  is  3.0. 

19.  The  trapezoidal  response  curve  from  figure  13  with  the  steep  fall  off  region  enlarged 
with  the  bias  point  at  which  the  calculation  was  performed. 

20.  The  distortion  calculations  showing  the  harmonics  for  a  single  tone  at  10  GFlz  bias 
at  midpoint  in  figure  19. 

2 1 .  The  distortion  calculations  showing  the  fundamental  and  the  third  harmonic  generated 
intermodulation  distortion  (IMD3),  with  the  different  values  marked  by  the  bias  points 
in  figure  19. 

22.  The  structure  of  the  polymer  waveguide  for  the  Mach-Zehnder  modulator,  with  the 
electrode  structure  modified  to  obtain  50  Q  impedance,  velocity  matched,  and  with 
10  pm  electrode  thickness. 

23.  The  variable  coupling  may  be  implemented  by  variable  spacing,  or  by  a  third  guide 
between  the  two  guides  of  the  coupler,  or  alternatively  by  varying  the  etching  between 
the  guides. 

24.  The  variable  coupling  may  be  implemented  by  variable  spacing  between  the  guides, 

and  for  the  basic  coupler  in  the  polymer  guides  with  a  spacing  of  10  pm,  and  the  range  of 
coupling  is  also  shown. 

25.  The  variable  coupling  may  be  implemented  by  a  third  guide  between  the  guides  of 
the  basic  polymer  coupler  with  a  spacing  of  10  pm,  and  the  range  of  coupling  is  also 
shown. 


2 

4 

4 


5 

6 
7 
7 


9 

9 

9 

10 

10 

11 

11 

12 

12 

13 

13 

14 

14 

15 

16 

16 


11 


26.  The  variable  eoupling  may  be  implemented  by  varying  the  etching  depth  between 
the  guides  of  the  basic  polymer  guide  coupler  with  a  spacing  of  10  pm,  the  range 

of  coupling  is  shown.  16 

27.  The  electrode  structure  with  the  even  microstrip  mode  for  the  polymer  coupler  with  its 

characteristic  impedances  shown.  17 

28.  The  electrode  structure  with  slow  wave  implementation  for  the  coupler  in  even  mode.  1 8 

29.  Microwave  index  for  electrode  structure  in  figure  28,  even  mode  excitation.  18 


List  of  Tables 


I  Parameters  for  different  materials  for  switehing  voltage  length  ealeulation.  5 

II  Voltage-length  produet  values  for  different  designs,  with  p  values  and  1/Le  6 

for  eoupler  designs. 


IV 


1.  Objective 

The  objective  of  this  study  was  the  linearization  of  optical  modulator  response  in  the  variable  coupling 
directional  coupler  modulator  with  particular  emphasis  on  the  realization  in  polymer  waveguide  structures. 


2.  Introduction 

The  co-directional  coupler  electro-optic  modulator  has  an  additional  degree  of  freedom  that  is  not  available 
with  the  traditional  Mach-Zehnder  interferometric  modulator.  This  additional  parameter  for  the  coupler 
modulator  is  the  coupling  function  which  can  vary  over  the  length  of  the  device,  and  the  variable  coupling 
form  may  be  synthesized  to  obtain  the  required  response  function.  Thus,  a  linear  response  modulator  may 
have  its  coupling  function  synthesized,  and  in  fact  a  device  of  this  type  has  been  realized  in  GaAs/AlGaAs 
[1]  previously.  The  synthesis  method  in  this  project  uses  the  inverse  Fourier  transform  method  as  a  first 
guess,  as  it  is  can  be  shown  that  the  coupling  function  and  the  output  guide  signal  field  are  approximate 
Fourier  transform  pairs  [2].  The  synthesis  then  is  performed  using  a  modified  form  of  Newton's  method 
with  this  first  guess  as  its  initial  value  [3].  In  this  project  the  technique  has  obtained  two  linear  modulator 
designs,  and  the  switching  voltages  of  these  are  compared  to  the  traditional  Mach-Zehnder  and  the  uniform 
directional  coupler  designs. 

The  coupling  functions  synthesized  by  this  technique  have  real  values  obtained  by  choosing  symmetric 
response  functions.  Flowever,  these  functions  also  have  negative  coupling  regions,  and  the  negative 
coupling  regions  are  obtained  by  placing  a  n  phase  shift  in  one  of  the  guides  relative  to  the  other. 
Subsequent  return  to  positive  coupling  is  obtained  by  inserting  another  phase  shift  region.  In  the  linear 
modulator  [2]  this  was  realized  by  increasing  the  length  of  the  one  of  the  guides  relative  to  the  other.  In  the 
case  of  the  GaAs/AlGaAs  structure  the  effective  index  is  about  3.4  at  1300  nm  wavelength,  and  the 
required  increase  is  0. 19  pm  in  one  of  the  guides.  In  the  case  of  the  polymer  guide  with  an  effective  index 
of  1.61,  this  increase  in  length  is  about  0.4  pm.  Flowever,  for  polymer  guides,  which  are  low  contrast 
structures,  the  bend  radii  for  insertion  of  this  extra  length  is  large,  of  the  order  of  millimeters  for  low  loss 
and  in  these  structures  a  different  approach  was  required.  The  technique  suggested  is  a  section  of  the 
coupler  in  which  the  two  guides  have  slightly  different  effective  indices,  and  by  choosing  an  appropriate 
length  the  phase  shift  of  n  may  be  readily  obtained. 

The  electrode  structure  of  the  current  polymer  guide  Mach-Zehnder  modulators  are  microstrip  lines  with 
the  Au  metallization  of  about  1  pm  thickness,  and  8  pm  wide.  The  skin  depth  in  Au  metallization  at  1  GHz 
is  about  3  pm,  which  leads  to  a  lossy  electrode  structure.  In  this  report  a  low  loss  Au  metallization  design 
for  the  existing  Mach-Zehnder  modulators  has  been  obtained.  For  the  coupler  modulator,  the  microstrip 
lines  on  the  two  guides  become  coupled  and  thus  even  and  odd  modes  may  be  obtained  by  the  form  of 
excitation.  In  this  case,  it  is  suggested  that  the  even  mode  be  used,  with  the  polymer  poled  in  opposite 
directions  under  the  two  electrodes  to  provide  the  necessary  push  pull  effect.  However,  in  this  case  the 
impedance  of  the  even  mode  microstrip  line  is  of  the  order  of  35  ohms,  including  the  slow  wave  structure 
for  matching  the  phase  velocity  and  the  optical  guide  velocity.  For  even  mode  excitation,  the  domains 
under  the  two  optical  waveguides  would  have  to  be  reverse  poled.  Using  the  odd  mode  excitation  for  this 
line  results  in  higher  impedance,  but  this  was  not  examined  due  to  time  constraints. 

The  following  sections  summarize  the  work  in  these  areas. 


3.  Design  Technique  for  Coupler  Modulators 

The  schematic  diagram  for  the  coupler  is  shown  in  Figure  1,  with  the  two  guides  in  close  proximity  to  each 
other,  with  the  input  to  guide  one  given  as  R(0)  and  the  input  to  the  second  guide  S(0).  When  a  signal  is 
placed  on  the  first  guide  of  the  uniform  coupler  modulator,  both  the  even  and  odd  supermodes  of  this 
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coupler  are  excited,  with  the  modes  adding  constructively  in  this  first  guide  and  destructively  in  the  second 
guide.  The  signal  which  is  now  the  superposition  of  the  two  modes  travels  down  the  coupler  with  the  odd 
mode  velocity  being  slightly  higher.  A  specific  distance  later  the  phases  of  the  odd  and  even  modes  are 
interchanged  because  of  the  difference  in  velocities,  so  that  the  modes  add  constructively  in  the  second 
guide,  but  destructively  in  the  first,  in  effect  causing  a  transfer  of  power  from  guide  one  to  guide  two. 


^R(L)^S(L) 


R(0) 

S(0) 

Figure  1 :  Schematic  of  the  uniform  coupling  coupler  modulator 


The  coupled  mode  equations 


for  the  directional  coupler  are  given  by: 

—  =  jdR-jK{z)Se-^*^^^ 
oi 

dz. 


(1) 


where  R  and  S  are  the  complex  amplitude  fields  in  the  guides  as  a  function  of  z,  spatial  position  along  the 
coupler,  k(z)  is  the  coupling  coefficient  which  is  a  slowly  varying  function  of  z,  Ap  is  the  difference  in  the 
propagation  constants  of  the  two  guides  uncoupled,  5  =  Ap/2,  and 

^(z)=  J[Ay0(z')-AAO)]*'  (2) 

0 

In  the  directional  coupler  modulator,  the  guides  are  made  identical  to  each  other,  so  that  5  is  zero  in  the 
unbiased  condition,  and  the  application  of  bias  makes  Ap  nonzero,  and  hence  5  is  also  nonzero. 


The  solution  to  these  equations  for  a  uniform  coupler  of  length  L,  with  initial  conditions  of  R{Q)  =\  and 
R(0)  =0,  is  given  by: 


i?(A)  =  cosfz/c-^l  +  {SjKf 


S{L) 


-ysin( 

Lk^ 

'iTtovfj 

A 

/l  +  (<5Al 

jS  sin[ 

Lk^i 

A 

/i + [sikY 

(3) 

(4) 


The  amplitude  response  of  the  uniform  coupler  modulator  is  given  by: 
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(5) 


77^S(L)\ 


2  _ 

R(0)=l,S(0)=0~ 


-j  sin^  (Lfc-\Jl  +  (S/fc)^ 
l  +  (S//cf 


Define  the  eoupling  length  Lc  as  the  distanee  where  all  the  power  in  the  first  guide  of  the  uniform  eoupler  is 
transferred  to  the  seeond,  or  |S|  =  1.  In  the  ease  when  there  is  no  bias  whieh  implies  d=  0,  then  from  the 
above  equation  L  =  L^  =  7I:I(2k).  For  the  ease  when  bias  is  applied  so  that  6  is  not  zero,  then  for 


L  =  4,  for  5  =  0,  Ik=^  or  =  ^[3 71:12.  Also  when  the  two  guides  are  biased  so  as 

to  have  different  effeetive  indiees,  then  by  definition 

^pl2  =  {n  -  n  )7tl2^  =  2An 


guide^^  ^0 


(6) 


The  differenee  in  the  guide  effeetive  indiees  due  to  bias  is  given  by  [2] 


An 


guide 


=  (n^rj^)l2d 


(7) 


where  n  is  the  guide  mean  effeetive  index,  r^x  is  the  eleetro-optie  eoeffieient  depending  on  the  direetion  of 
the  field  interaeting  with  the  optieal  mode,  V  is  the  voltage  applied  to  the  eoupler  eleetrodes,  d  is  the 
distanee  between  the  eleetrodes,  and  F  is  the  overlap  integral  between  the  bias  eleetrie  field  and  the 
eleetrie  field  of  the  optieal  mode.  Substituting  for  the  $\delta$  from  the  above  equation  : 


VswirchLc=pUQd)/(n^rT)  (8) 

where 

F  =  =  (9) 


for  push  pull  eleetrode  exeitation  of  the  uniform  eoupler. 


Using  the  transformation  to  the  Rieatti  equation  [1,2],  it  ean  be  shown  for  very  small  eoupling  that  Sand  K 
are  related  through  the  Fourier  transform  equation  of  the  form: 

k(z)  =  —  }  I  S(^  I  cos,{2dz)dd  (10) 

^  0 

Thus,  if  the  desired  intensity  response  for  the  eoupler  modulator  is  given  by  |  S(S)  f ,  the  first  guess  with 
small  eoupling  is  obtained  from  the  above  equation  (10).  The  eoupling  funetion  derived  from  this  equation 
extends  to  ±go  ,  and  needs  to  be  truneated  symmetrieally  to  prevent  ehirp  being  generated.  The  tmneated 
eoupling  funetion  will  only  provide  an  approximate  form  of  the  desired  response,  and  thus  needs  to  be 
modified  to  provide  a  eloser  agreement.  The  next  step  uses  the  Newton's  method  [3],  with  this  truneated 
response  as  the  first  guess.  The  error  in  the  response  funetion  is  evaluated,  the  appropriate  Jaeobian  is 
obtained,  and  by  the  least  squares  method  the  eorreetion  is  obtained.  Subsequent  iterations  move  the 
eoupling  funetion  eloser  to  the  desired  response  funetion.  Figure  2  shows  the  desired  response  funetion  in 
the  form  of  a  trapezoid  and  the  iterated  response  from  the  modified  eoupling  funetion  is  shown  in  Figure  3. 


3 


Figure  2:  The  desired  trapezoidal  response  function  obtained  after  several  iterations,  the  normalized 
intensity  response  plotted  against  the  normalized  switching  voltage. 


Figure  3:  The  iterated  coupling  function  for  the  response  function  in  figure  2. 

The  power  in  the  second  guide  along  the  coupler  is  shown  in  Figure  4.  Note  that,  for  the  zero  bias  case, 
shown  in  blue,  the  transfer  of  power  to  the  second  guide  only  takes  place  in  the  major  lobe  of  the  coupling 
function  in  the  center  of  the  coupler.  Thus,  the  regions  of  small  coupling  at  the  ends  of  the  coupling 
function  contribute  to  shaping  the  response  function.  The  green  curve  shows  the  power  in  the  guide  with 
bias  halfway  down  the  trapezoidal  slope  region,  ft  becomes  apparent  that  the  power  is  initially  transferred 
to  the  second  guide,  and  then  half  of  it  returns  to  the  first  guide  to  obtain  a  final  value  of  half  The  red 
curve  shows  the  modulator  biased  to  zero  output  at  the  second  guide,  with  initial  transfer,  and  then  the 
power  returning  entirely  to  the  first  guide. 
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Figure  4:  Plot  of  normalized  power  against  normalized  distance  along  the  coupler  for  different  bias  values. 
The  blue  curve  is  for  zero  bias,  green  curve  is  for  bias  which  obtains  half  power,  and  the  red  curve  is  for 
bias  when  no  transfer  takes  place  into  the  second  guide. 


4.  Calculation  of  Switching  Voltage-Length  Product  for  Different  Modulator 
Designs 


When  the  guide  electrodes  of  a  modulator  are  biased,  the  index  in  the  guides  vary  resulting  in 
becoming  non-zero,  and  in  turn  this  results  in  d  becoming  nonzero.  As  shown  above  for  the  push  pull 
bias,  the  product  of  the  switching  voltage  V switch  the  modulator  length  is  given  by  the  equation: 


^ switch^  c  P 


AqcI 

n^rT 


(11) 

where  p  is  given  by  /  TT ,  and  is  the  figure  of  merit  for  the  designs,  F  is  the  overlap  integral 

between  the  optical  field  and  the  bias  field,  d  is  the  electrode  gap,  n  is  the  effective  index,  is  the 
electro-optic  coefficient  for  the  geometry  being  used.  For  the  Mach-Zehnder  interferometric  modulator  p  is 
0.5,  and  for  the  constant  coupling  coupler  modulator  p  is  0.87.  Table  I  provides  the  electro-optic  constants 
for  different  materials  used  for  the  switching  voltage-length  calculations,  together  with  other  parameters. 


Table  I:  Parameters  for  different  materials  for  the  switching  voltage-length  calculation. 


Parameters 

CLD-75  Polymer 

G3A.s/Alo.2Gr3o,sAs 

LiNiBOj 

Wavelength  (nm) 

1550  nm 

1550  nm 

1550  nm 

Electrode  Spacing  d  (pm) 

8  pm 

2  pm 

8  pm 

Overlap  integral  F 

0.95 

0.65 

0.55 

Electro-optic  coefficient  rxx 

36pm/V 

1.43  pm/V 

31.5  pmW 

Refractive  index  n  at  1550 

nm 

1.612 

3.43 

2.138 

n  I'xx 

150 

58 

308 

Ratio  (VswitchL)/  p 

8.66  V-cm 

8.26  V-cm 

7.32  V-cm 
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Table  II:  Voltage-length  produet  values  for  different  designs,  with  p  values,  and  L/Lc  for  eoupler  designs. 


Modulator  Design 

P 

L/Lc 

CLD-75 

Polymer 

GaAs/ Alo.2Gao.8As 

LiNiBOs 

Uniform  coupler 

0.87 

1 

7.5  V-cm 

7.15  v-cm 

6.3  V-cm 

Mach-Zehnder 

0.5 

4.33  V-cm 

4.13  V-cm 

3.66  V-cm 

Coupler:  Trapezoidal 
response 

1.2 

1 

10.4  V-cm 

9.9  V-cm 

8.8  V-cm 

Coupler:  Linear  response 

4.1 

1 

35.5  V-cm 

33.9  V-cm 

30.0  V-cm 

Coupler:  constant  k,  n  phase 
shifts 

1.8 

2 

15.6  V-cm 

14.9  V-cm 

13.2  V-cm 

Coupler:  5  varied,  constant  k 

2.7 

3 

23.4  V-cm 

22.3  V-cm 

19.8  V-cm 

The  spaeing  for  the  GaAs/Alo  2Gao.8As  modulator  with  Sehottky  eontaets  has  all  the  voltage  drop  in  the 
depletion  layer.  Although  the  eleetrode  spaeing  may  be  8  pm  to  be  eonsistent  with  the  other  eases,  a 
realistie  value  here  depends  on  the  intrinsie  unintentional  doping  level.  This  was  ehosen  to  be  2  pm  on 
average  for  both  MBE  and  MOCVD  material.  Note  that  the  semieonduetor  estimate  is  at  best  approximate 
sinee  the  depletion  layer  thiekness  varies  as  the  square  root  of  the  sum  of  the  bias  voltage  and  the  Sehottky 
barrier  built-in  voltage. 

Congruent  LiNiB03,  as  opposed  to  the  stoichiometrie  LiNiB03,  is  widely  used  and  its  value  of  r3i  was  used 
in  the  above  Table  I.  The  overlap  integral  F  has  been  widely  quoted  as  0.20  for  lithium  niobate,  and  our 
earlier  estimate  of  the  ratio  (Vswitch  L)/  p  in  this  material  obtained  a  value  of  17.9  V-em  with  this  value  of 
r  .  However,  recent  discussions  with  Dr.  SriRam  of  Srico  Inc.  has  provided  us  with  a  better  estimate  of 
0.55  for  the  z-cut  structures  and  this  has  resulted  in  the  ratio  (Vswitch  L)/  p  shown  in  Table  I. 

The  various  designs  that  have  been  examined  in  this  study  are  listed  in  Table  II  together  with  their  Vswitch  L 
values  in  the  appropriate  material  system.  In  the  following  figures  the  different  designs  are  shown  with  the 
corresponding  intensity  versus  bias  voltage  response  curves.  Figures  5  and  6  show  the  schematic  diagram 
of  the  Mach-Zehnder  interferometric  modulator  and  its  calculated  response  function  with  a  straight  line 
linear  response  showing  the  deviation  respectively.  The  calculated  value  of  p  is  0.5  for  the  push  pull  case. 

Figure  7  shows  the  schematic  diagram  of  the  uniform  coupler,  and  Figure  8  shows  the  response  function, 

the  value  of  p  for  this  device  is  0.866  /2 ) .  The  linear  response  coupler  case  obtained  the  coupling 

from  the  Fourier  transform,  as  shown  in  Figure  9,  and  the  corresponding  response  in  Figure  10  for  the 
truncated  coupling  function. 

The  trapezoidal  response  coupler  design  coupling  function  is  shown  in  Figure  1 1  and  the  response  function 
in  Figure  12,  desired  and  calculated.  Since  the  realization  of  this  coupling  function  is  difficult,  the 
alternative  coupling  function  with  flat  topped  regions  and  the  corresponding  response  function,  which  has 
suffered  from  this  modification,  are  shown  in  Figures  13  and  14. 


-b  V 


-V 

Figure  5:  Schematic  diagram  of  the  Mach  Zehnder  modulator 
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Figure  6:  The  Mach  Zehnder  calculated  response. 


-V 

Figure  7 :  Schematic  diagram  of  the  unifomi  coupler. 
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Figure  8:  The  response  of  the  uniform  coupler  in  which  the  normalized  transmitted  power  is  plotted  against 
the  normalized  voltage.  The  value  of  p  for  this  coupler  is  VS  / 2  for  the  push  pull  case. 
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Figure  9:  The  linear  response  coupler  modulator  coupling  function  obtained  from  the  Fourier  transform  of 

the  desired  linear  response. 
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Figure  10:  The  linear  response  coupler  modulator  showing  the  response  function  as  obtained  from  the 

coupling  function  in  figure  9. 


Figure  1 1 :  The  trapezoidal  response  coupler  design  showing  the  coupling  function. 


Figure  12:  The  trapezoidal  response  coupler  design  showing  the  response  function, 

desired  and  calculated. 
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Coupling  Function  for  Trapezoid  Response 


Figure  13:  The  trapezoidal  response  eoupler  design  from  figure  11  showing  the  modified  eoupling  function 
_ to  have  flat  top  regions. _ 


Figure  14:  The  trapezoidal  response  coupler  design  from  figure  13  showing  the  resultant  response  function. 


Two  other  designs  are  also  considered.  In  the  first  design,  the  coupling  is  made  positive  or  negative  but 
retains  a  constant  magnitude.  This  is  to  be  achieved  by  placing  1 80  degree  phase  shifts  at  the  sign  change 
positions,  and  this  is  discussed  further  in  the  following  section.  The  second  design  maintains  k  constant  but 
switches  the  sign  of  5  in  the  biased  condition,  either  by  having  segmented  electrodes  which  are  biased 
alternately  with  opposite  polarities  or  by  poling  the  material,  polymer  or  LiNiB03  or  other  ferroelectric 
material  in  alternate  directions.  In  the  former,  the  effect  of  the  sign  change  for  k  results  in  a  slight  dip  in 
the  transfer  of  power  along  the  coupler,  but  ultimately  results  in  a  very  linear  response  function  as  shown  in 
Figure  15  and  Figure  16  for  a  longer  device  of  twice  the  coupling  length,  with  p=1.8.  In  the  latter  case,  the 
change  in  bias  is  similar  to  the  delta-beta  coupler  modulators,  but  also  allows  the  use  of  poling  to  obtain 
this  voltage  bias  change  with  continuous  electrodes.  In  this  case  adjusting  the  positions  and  the  lengths  of 
the  poled  sections  results  in  a  linear  response  modulator  as  shown  in  Figures  17  and  18  with  the  device 
length  three  times  the  coupling  length,  and  p  having  a  value  of  2.7. 
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Coupling  Function  for  Linearized  Response 
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Figure  15:  The  constant  |k|  coupler  modulator  design  but  with  sign  changes  for  the  coupling. 


Figure  16:  The  constant  |k|  coupler  modulator  but  with  sign  changes  for  the  coupling.  The  response 
function  is  shown,  the  value  ofp  is  1.8  and  the  L/Lc  ratio  is  2.0. 
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Figure  17:  The  eonstant  |k|  coupler  modulator  design  but  with  sign  changes  for  5  along  the  length  of  the 

coupler. 


Low-Voltage  Response  with  Varied  Delta  Polarity 


Figure  18:  The  response  function  of  the  constant  |k|  coupler  modulator  but  with  sign  changes  for  5  along 
the  length  of  the  coupler,  the  value  of  p  is  2.7  and  L/Lc  is  3.0. 


5.  Distortion  Calculations 


The  usual  technique  for  calculating  the  distortions  in  modulators  is  to  express  the  response  function  as  a 
polynomial  in  the  applied  voltage,  and  thus  for  the  present  case  the  output  signal  power  is  expressed  as: 

1^12 


\out-  ^0  +  ^y^ppl  +  ^^appl 


+  . 


(12) 

Flowever,  expressing  the  response  function  in  polynomial  form  results  in  considerable  error  and  the  results 
are  not  accurate. 


Another  method  is  to  superpose  two  signals  at  slightly  different  frequencies,  sum  them  and  take  the  Fourier 
transform,  either  in  the  FFT  technique  or  in  the  usual  integral  form.  Flowever,  it  is  necessary  that 
periodicity  is  not  violated  by  taking  too  small  a  sequence,  and  this  is  also  results  in  inaccuracies  due  to  the 
large  number  of  points  required.  The  most  accurate  method  is  to  use  the  2D  FFT  technique  in  which  the 
two  dimensions  correspond  to  the  time  periods  of  the  two  tones  passed  through  the  modulator 
simultaneously.  Matlab  and  other  programs  have  this  as  a  standard  function,  but  it  is  also  possible  to 
program  it  to  work  with  the  single  FFT  code.  Figure  19  shows  the  trapezoidal  response  function,  with  the 
enlarged  section  of  the  response  curve.  Figures  20  and  21  show  the  harmonics  and  the  intermodulation 
distortion  third  harmonic  (IMD3),  which  are  2fi-  f2  and  2f2-  fi  components  with  different  drive  levels. 
Drawing  the  noise  floor  provides  the  spur  free  dynamic  range  (SFDR),  which  is  the  vertical  distance  from 
the  intersection  of  the  noise  floor  of  the  IMD3  line  to  the  fundamental  power  at  that  point. 
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T  rapezoidal  R  espon  se 


Normalized  Bias 


Figure  19:  The  trapezoidal  response  eurve  from  figure  13  with  the  steep  fall  off  region  enlarged  with  the 
bias  point  at  whieh  the  calculation  was  performed. 


Figure  20:  The  distortion  calculations  showing  the  harmonics  for  a  single  tone  at  10  GFlz  bias  at  midpoint 
in  figure  19. 


13 


Fiiii(l<inieiiti')l  iIikI  Hiiniionics 


Figure  21:  The  distortion  calculations  showing  the  fundamental  and  the  third  harmonic  generated 
intermodulation  distortion  (IMD3),  with  the  different  values  marked  by  the  bias  points  in  figure  19. 
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Figure  22:  The  structure  of  the  polymer  waveguide  for  the  Mach-Zehnder  modulator,  with  the  electrode 
stmcture  modified  to  obtain  50  Q  impedance,  velocity  matched,  and  with  10  pm  electrode  thickness. 

6.  Linearized  Coupler  Designs  in  Polymer  Waveguides 

The  polymer  modulator  guide  structure  is  shown  in  Figure  22  from  the  work  of  the  USC-UCLA  groups  [4], 
except  that  the  width  of  the  electrode  structure  is  12  pm  and  the  thickness  is  10  pm,  as  per  calculations 
performed  in  this  project  for  velocity  matched  50  Q  microstrip  line.  The  original  design  for  their  Mach- 
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Zehnder  modulator  had  8  |im  wide  electrodes,  1  |im  thick,  leading  to  a  velocity  matched  microstrip  line 
with  impedance  of  70  Q  for  the  individual  lines.  The  polymer  guides  are  poled  with  bias  on  the  electrodes 
in  the  Mach-Zehnder  modulator,  and  in  the  coupler  it  is  suggested  that  this  type  of  poling  is  performed  in 
which  the  two  guides  are  in  oppositely  poled  material.  Since  our  design  suggest  a  10  pm  spacing,  it  will  be 
necessary  to  ensure  that  surface  breakdown  does  not  prevent  poling  to  take  place.  If  this  cannot  be 
performed  successfully,  then  we  will  have  to  alter  our  electrode  design  from  the  coupled  even  mode  slow 
wave  structure  to  the  odd  mode  slow  wave  structure. 


6.1  Realization  of  Variable  Coupling  Designs  in  Polymer  Waveguides 

6.1.1  Variable  Coupling  Schemes 

The  variable  coupling  designs  may  be  implemented  using  alternative  approaches  as  shown  in  Figure23. 


Figure  23:  The  variable  coupling  may  be  implemented  by  variable  spacing,  or  by  a  third  guide  between  the 
two  guides  of  the  coupler,  or  alternatively  by  varying  the  etching  between  the  guides. 


The  simplest  method  is  the  variable  spacing  scheme,  however,  the  problem  of  electrodes  following  the 
variable  spacing  remains,  and  since  this  changes  the  electrode  impedance  this  needs  to  be  accounted  for. 
Beam  propagation  code  from  RSoft  was  used  to  determine  the  range  of  coupling  that  may  be  implemented 
with  this  method,  and  the  results  are  shown  in  Figure  24.  These  results  show  that  this  is  the  simplest 
method  of  implementing  the  variable  coupling  concept,  with  a  wide  range  of  coupling  that  may  be 
implemented,  and  a  smoothly  varying  coupling  function  may  also  be  implemented. 


The  second  scheme  is  to  place  a  third  guide  between  the  guides  of  the  coupler.  In  this  case,  the  third  guide 
needs  to  be  single  mode  to  ensure  that  no  multimode  propagation  occurs,  and  also  it  is  necessary  to  retain 
the  basic  odd  and  even  mode  structure  of  the  coupler,  so  that  only  this  basic  eigenvector  of  this  three  guide 
coupler  is  excited.  The  range  of  coupling  that  may  be  obtained  with  this  method  is  liiuited  to  a  two  to  one 
ratio  as  shown  in  Figure  25,  obtained  with  the  RSoft  BeamProp  program.  The  calculations  show  that  a 
smoothly  varying  coupling  function  may  also  be  implemented  within  this  range  of  variation. 


The  third  scheme  is  to  have  the  gap  between  the  coupled  guides  constant,  but  to  vary  the  etch  depth 
between  them.  While  a  wide  range  of  coupling  may  be  obtained  with  this  method,  the  problem  is  that  when 
the  cladding  surface  approaches  the  guide  layer,  the  mode  becomes  very  lossy,  due  to  the  etch  roughness. 
Figure  26  shows  this  for  the  polymer  coupler,  again  using  BeamProp  for  the  calculations. 
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Computed  Transverse  Mode  Profile  (m=0,n^jj=1. 597021) 
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Figure  24:  The  variable  coupling  may  be  implemented  by  variable  spacing  between  the  guides,  and  for  the 
basic  coupler  in  the  polymer  guides  with  a  spacing  of  10  pm,  and  the  range  of  coupling  is  also  shown. 


Computed  Transverse  Mode  Profile  (m=0,n  „=1 .597054) 
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Figure  25:  The  variable  coupling  may  be  implemented  by  a  third  guide  between  the  guides  of  the  basic 
polymer  coupler  with  a  spacing  of  10  pm,  and  the  range  of  coupling  is  also  shown. 


Dual  Polymer  Waveguide  Coupling  Length  vs.  Etch  Depth 


Figure  26:  The  variable  coupling  may  be  implemented  by  varying  the  etching  depth  between  the  guides  of 
the  basic  polymer  guide  coupler  with  a  spacing  of  10  pm,  the  range  of  coupling  is  shown. 


16 


6.1.2  Negative  Coupling  by  Phase  Shifts 

Negative  eoupling  may  be  implemented  by  plaeing  a  relative  phase  shift  between  waveguides.  This 
requires  one  waveguide  to  be  A,  12  longer  than  the  other.  In  semieonduetor  waveguides,  in  our  previous 
work  [3],  the  gap  was  inereased  between  the  guides  and  then  an  inereased  A 12  length  was  added  to  one  of 
the  guides,  in  the  form  of  a  seetion  of  a  loop,  and  the  loop  radii  were  in  the  region  of  100  pm  or  larger.  In 
polymer  guides,  this  requires  bends  with  0.48  mm  or  larger  radii,  whieh  makes  these  rather  long.  The 
alternative  is  to  use  a  polymer  eoupler  with  different  waveguide  refraetive  indiees  to  produee  this  phase 
shift.  Thus,  leaving  one  guide  unpoled  produees  ISOdegrees  phase  shift  in  40  pm.  Implementation  of  this 
is  thus  very  simple  and  does  not  require  the  loop  phase  shifter. 


6.1.3  Electrode  Structure  for  Polymer  Coupler 

The  skin  effeet  thiekness  of  gold  eleetrodes,  traditionally  used  in  modulators  at  1  GHz  is  about  3  pm  and 
therefore  it  is  neeessary  to  have  eleetrodes  with  thiekness  of  10  pm  to  have  minimum  eonduetor  loss  for 
modulation  at  10  GHz  upwards.  The  polymer  guides  use  mierostrip  lines  for  their  eleetrodes  and  the 
widths  are  adjusted  so  that  the  phase  veloeity  is  elose  to  the  veloeity  of  light  in  the  guides,  e/1.6.  In  the 
ease  of  the  polymer  eoupler,  eoupled  mierostrip  lines  are  used  in  the  even  mode,  so  that  the  optieal  guide 
regions  under  the  two  guides  are  poled  in  opposite  direetions.  With  10  pm  thiek  eleetrodes  eonsiderable 
energy  propagates  in  the  air,  and  so  these  beeome  fast  wave  struetures,  and  to  veloeity  mateh  them  we  need 
to  introduee  additional  eleetrodes  whieh  form  slow  wave  struetures.  However,  the  impedanee  suffers,  as 
this  is  now  in  the  35  Q  range.  These  are  shown  in  Figures  27,  28  and  29. 
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Figure  27:  The  eleetrode  stmeture  with  the  even  mierostrip  mode  for  the  polymer  eoupler  with  its 

eharaeteristie  impedanees  shown. 
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Figure  28:  The  electrode  structure  with  slow  wave  implementation  for  the  coupler  in  even  mode. 
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Figure  29:  Microwave  index  for  electrode  structure  in  figure  28,  even  mode  excitation. 


7.  Conclusions 

The  project  conducted  a  study  of  two  aspects  of  the  variable-coupling  linear  response  directional  coupler 
modulator.  The  first  part  of  the  study  detemiined  the  optimum  method  of  synthesizing  the  coupling 
function  for  linear  response  in  conjunction  with  the  implementation  constraints.  The  second  part  of  the 
study  examined  the  problems  involved  in  implementing  the  variable-coupling  linear  response  directional 
coupler  modulator  in  polymer  waveguides.  In  conclusion,  this  study  shows  that  variable  coupling 
directional  coupler  modulators  with  ultra-linear  response  may  be  readily  realized  in  polymer  waveguides. 
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